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Abstract 

The selection rules for dipole and Raman activity can be relaxed 
due to local distortion of a crystalline structure. In this situation 
a dipole-inactive mode can become simultaneously active in Raman 
scattering and in dipole interaction with the electromagnetic field. 
The later interaction results in disorder-induced polaritons, which 
could be observed in first-order Raman spectra. We calculate scatter- 
ing cross-section in the case of a material with a diamond-like average 
structure, and show that there exists a strong possibility of observing 
the disorder induced polaritons. 

PACS numbers: 78.30-j, 63.50+x, 63.20-e. 

1 Introduction 



The effects of a resonance linear interaction between excitations of differ- 
ent nature are well known and their importance for properties of crystalline 
materials has been appreciated for a long time. One could mention, for in- 
stance, such phenomena as magneto-elastic resonance QXQ , or various kinds of 
the polariton resonance |2|. The former refers to resonance interaction be- 
tween spin and elastic waves that results in mixed excitations carrying both 
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magnetic and elastic types of energy. The latter term is used to refer to exci- 
tations arising as a result of interaction between electromagnetic waves and 
excitations of crystals such as phonons, excitons, plasmons, and so on. In 
spite of significant differences between all the cited examples, they share one 
common feature. In the absence of interaction the dispersion curves of the 
participating excitations cross each other at a certain point. This crossing is 
actually responsible for the resonance nature of the interaction, which lifts 
the degeneracy of the initial excitations at the crossing point and gives rise 
to the mixed states. 

In a recent paper [Q] the effect of crossing resonance has been studied 
in systems with a random interaction between participating excitations. It 
was shown that even a random coupling parameter with zero mean value 
can, under certain conditions, result in the strong coupling between interact- 
ing waves, leading to a lifting of the degeneracy at the crossing point. The 
mixed excitations arising in this model of disorder-induced crossing resonance 
(DICR), however, differ significantly from the respective excitations in the 
deterministic case. First of all, excitations that are actually mixed together 
are coherent excitations of one nature with nonzero average amplitude and 
non-coherent scattered excitations of the other nature with average amplitude 
equal to zero. Correspondingly, one has to consider two such mixed excita- 
tions, and each of them can result in dispersion laws split at the crossing point 
in two branches . It is possible, therefore, to have three or even four differ- 
ent dispersion curves describing averaged properties of the system. Recently 
the model of DICR was developed further to include nonzero mean value of 
the random coupling parameter J|. It was shown in Ref.[|J that transition 
between the completely disordered case with zero mean to the nonrandom 
situation occurs in a nontrivial way. The idea of the DICR was first applied 
to the magnetoelastic resonance in amorphous zero-mean-magnetostriction 
alloys || . In magnetoelastic case it is possible to perform a direct observation 
of elastic and magnetic susceptibilities modified due to DICR. The relation- 
ship between maxima of the susceptibilities and solutions of the dispersion 
equations for magnetoelastic resonance was studied in Ref. |]] 

In Ref.0 it was suggested that DICR can also occur for the polariton 
resonance in disordered materials if one considers the interaction of electro- 
magnetic waves with phonon or exciton excitations that would be dipole- 
inactive in the absence of disorder (disorder- induced polaritons). Polariton 
dispersion curves in ideal crystals can be experimentally studied by means of 
Raman scattering (see, for example, Ref. 0). It was also suggested in Ref.0 
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that polaritons arising due to retardation effects in the electromagnetic in- 
teraction between local impurity vibrations can also be observed in Raman 
scattering. It seems natural, therefore, to explore whether it is possible to 
use the Raman effect to observe the disorder-induced polaritons. In this pa- 
per we present calculations of the first order Raman scattering cross-section 
from forbidden polaritons arising in a simple cubic diamond-like structure 
due to local structural distortions. 

2 Intensity of Raman scattering from disorder- 
induced polaritons. General expression 

The polarization, P, of a dielectric medium subjected to an external elec- 
tromagnetic field, E in ,with frequency u in can be expressed by means of a 
susceptibility tensor Xij as 

Pi = XiMn- (1) 

In the frequency region of our interest, susceptibility can be divided into 
two major contributions: ion susceptibility Xij n an d electron susceptibility 
Xtj- The former part of the susceptibility arises due to the coherent mo- 
tion of ions caused by the incident electromagnetic wave and is responsible 
for a modification of the wave's velocity (real part) and for its absorption 
(imaginary part). The electron contribution xtj reflects the contribution of 
electron transitions in the dipole moment of the atoms constituting crys- 
tal. For frequencies far enough from those of electron transitions and in the 
linear approximation xtj can be considered as a frequency independent con- 
stant. However, if one takes into account nonlinear corrections such as an 
interaction between electrons and ionic vibrations, the electron contribution 
becomes modulated by thermal ionic vibrations giving rise to the inelastic 
scattering of light known as Raman scattering. As a result of this scattering 
the frequency, u in , and wave number k in , of the incident wave changes. In 
the case of scattering from an ideal crystal, frequencies and wave numbers of 
the incident and scattered waves obey well-known kinematic relations that 
reflect energy and momentum conservation laws: 

UJ in -U s = ±tt (2) 
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Here u s and k s are the frequency and wave number of the scattered field, and 
Q and K are the frequency and wave number of the excitations of the medium 
responsible for scattering. The signs + and — correspond to processes of 
emission and absorption of the excitations respectively giving rise to Stokes 
and anti-Stokes components of the scattered field. In the presence of disorder 
the second of the equations (0) does not hold any more since disordered 
systems lack translational invariance. 

Intensity of the field scattered at a given direction within a certain fre- 
quency interval is determined by the differential cross-section d 2 a/dodQ, 
where do is a differential of a solid angle at a given direction. In general 



the cross section can be presented as |10 
d 2 a 



dodVt 



const x e\(£efe? « 5 Xik (^, K)6 Xj i(£l, K) » (3) 



where e s and e m are polarization vectors of the scattered and the incident 
fields respectively, and the const in front of the whole expression is a combi- 
nation of such parameters as volume of the scattering region of the sample, 
frequency of the incident wave, etc., which are not significant for the pur- 
pose of our consideration. Double angular brackets << • • • >> denote two 
types of averaging: thermodynamical average over statistical ensemble, and 
the average over realizations of the random inhomogeneities of the system. 
Modulation of the susceptibility S Xi k(fi,K.) can be found if one expands 
the electronic susceptibility into a power series with respect to the ampli- 
tudes of excitations responsible for the scattering. In the case of scattering 
due to polaritons, the electronic susceptibility is modulated by both ionic 
displacements W a and electric field E associated with those displacements. 
Therefore, 5 X ik(Q, K) is the sum of two terms 

5 Xik (n, K) = £/ d 3 qW:(n, qR fc (K - q) + J d 3 qE* m (Q, q)^(K - q), 

(4) 

where W a represents an amplitude of cr-th normal mode of the ionic vibra- 
tions, h tk is the tensor of mechano-optic coefficients, and g tkm represents 
electro-optic coefficients of the system, wave number K is determined by the 
Eq.(||), and an additional integration over q reflects the inhomogeneous na- 
ture of the considered system and the lack of the translational invariance. 
For this reason susceptibility derivatives h tk and g lkm should be considered 
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as random functions. Within the phenomenological approach utilized in this 
paper we have to make additional assumptions about the properties of these 
functions based upon some general arguments. We begin with separating out 
mean parts of these functions and represent them in the following form: 



hf(K- q ) = <hf>5(K- q ) + hf(K- q ) (5) 
g ikm (K - q) = < g ikm > 8{K - q) + g ikm (K - q 



According to this representation the scattering cross-section has two different 
contributions: the coherent one, resulting from averaged parts of the coeffi- 
cients h and g, and the non coherent contribution due to inhomogeneous parts 
of the tensors. The main physical difference between these contributions is 
that the former conserves momentum (this is reflected by the presence of the 
5-function in the respective term of Eq.(|5])), and the latter does not. The 
properties of the mean parts of the susceptibility derivatives is determined 
by the average symmetry of the system. We assume the cubic average sym- 
metry of the crystal structure that preserves central symmetry, and that the 
normal vibrational mode of the system is a triply degenerate (in the absence 
of interaction with electromagnetic field) mode of even parity. (An example 
of such a system can be the mode in a diamond-like structure). This 
mode is dipole inactive on average and its coupling to the electromagnetic 
field occurs only due to random deviations of structure from the average one. 
This assumption ensures that we actually deal with the disorder-induced po- 
laritons, and at the same time it guarantees that < h %k > is not equal to 
zero since it is well known Jl0| that infra-red activity and Raman activity 
are mutually exclusive in the crystals with a center of inversion. The tensor 
structure of < h lk > in the case of a diamond like structure can be presented 



in the following form 1 1 



< hf >d™ = h\e lkm \ (6) 

where d J a is a polarization vector of a given mode, and \eikj \ is the Levi-Civita 
symbol. 

The question about properties of the electro-optical coefficient in our 
model is more subtle. In the regular case of crystalline materials these coef- 
ficients come into play only in the case of vibrational modes of odd parity, 
because only in this case does there exist an electromagnetic field associated 
with the vibrations. Then a third-rank tensor g tkm in such a situation is 
equal to zero in crystals with central symmetry, therefore, the whole issue of 
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the properties of the electro-optical coefficients in crystalline materials arises 
only in non-centro-symmetrical crystals. In our situation of disorder-induced 
polaritons the mode one deals with has even parity. Therefore, symmetry 
arguments do not forbid nonzero value of the tensor g tkm even in materials 
with central-symmetry average structure. Therefore we can assume for the 
tensor structure of g lkm the same form as in Eq. (|6|) . 

We can conclude, therefore, that the coherent contribution in the scat- 
tering cross-section consists of two terms 

= const x [n(0) + 1](J^ + Jg*), (7) 
where the vibrational contribution Iffi is 

It = h 2 \e lkm \\e jlp \ete]eTeTY,€dll^ [^(fl, K)] , (8) 



and the electromagnetic contribution I c £ h is given by the similar expres- 
sion 

It = g 2 h km \\e jlp \e^el n erim [^(fl, K)] . (9) 

Writing down these expressions we make use of the fluctuation-dissipation 
theorem and express the vibrational and electric field correlation functions 
in terms of thermodynamical Bose-factor n(Q) + 1, and the respective vibra- 
tional and electromagnetic linear-response functions rj^(Q, K) and t]f-(Q, K) 
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In order to describe non-coherent contributions in the cross-section we 
have to make assumptions regarding correlation properties of the suscepti- 
bility derivatives h l £ and g lkm . It is natural to assume that their fluctuations 
are both homogeneous and isotropic. These assumptions allow one to intro- 
duce the spectral densities (spatial Fourier-transforms) of the corresponding 
correlation functions and present them in the following form 

< ^(k 1 )(^(k 2 ))* >= ^'5(k 1 )5(k 1 - k 2 ) 

< ^(kOGT J (k 2 ))* >= GfJSikJS^ - k 2 ) 

where 5"(k) is a scalar spectral density of the structural inhomogeneities, 
if^and G^™ are isotropic tensors with a structure reflecting properties of 
coefficients h l £ and g lkm . 
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Calculations of the non-coherent contributions in the cross-section re- 
quires averaging of the expressions like 

« ^(kOW^ka^K - kx)(^(K - k 2 ))* » 

. Doing so we use the simplest decoupling of the correlators of this from as 
follows: 

< < W^W^hfiK - kx)(^(K - k 2 ))* »~ 

< < W^W;^) »< hf(K - kO(^(K - k 2 ))* > 

This approximation is justified in the case of sufficienlty weak inhomogeneities, 
which is assumed throughout the paper. Within this approximation the non- 
coherent cross-section can again be separated into two terms I^ oh and I™ coh 
with quantities I having the meaning similar to that in Eq.(|7|). It is natural 
to assume isotropy of the correlation tensors H i anc l G 1 ^ with additional 
symmetry properties following from the symmetries of the initial mechano- 
optical and electro-optical coefficients. For the vibrational contribution Iy^ oh 
in the non coherent scattering one then has 

Iff* = [ft + H 2 (e s e m ) 2 } £ J d 3 qlm [ V ^(Q, q)] S(K - q), (10) 



where Hi and if 2 are the only independent elements of the tensor H 3 . An 
expression for the non-coherent electromagnetic contribution is somewhat 

iikm 
T jnl 



more cumbersome because of more complicated tensor structure of Gi* 171 - 



jncoh 



[ Gl + G 2 (eV") 2 ] J d 3 qlm [ri(Q, q] S(K - q) + 
[Gs {ere? + e°e°) + G 4 (efef + efe*)] f d 3 qlm [^(O, q)] S(K - q) 



3 Linear-response functions for disorder-induced 
polaritons 

This section of the paper concerns with linear response functions rj^(Q, q) 
and T]ij(Q, q), which determine Raman scattering cross-section in our sys- 
tem. We start from dynamical equations describing ionic vibrations in the 
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following form 



+ S r i)^( r i) = ^E\v). (12) 

This equation must be complemented by Maxwell's equations 

H + ( ?VXVXE=^ (13, 
with polarization density P defined as 

P=^$> a U a + Xe «E. (14) 

Indices label ions within one elementary cell, U l is a Cartesian component 
of a displacement vector, D^Jr, r x ) is a dynamical matrix describing inter- 
action between ions of masses m a and charges e a . For given radius- vectors 
r, r 1 the dynamical matrix can be diagonalized with respect to sub-cell and 
Cartesian indices. If c l aa .(r) is a matrix that diagonalizes the dynamical ma- 
trix, the equation of motion ([12]) and the equation for polarization (|14]) can 
be rewritten in the following form 

f) 2 W 

+ J2^(r,r l )W a (r 1 ) = Z„(r)4#(r) (15) 

ri 

P=^^^(r)d (T W (7 + Xei E, 

where normal mode W a {r) and its effective charge Z a (r) are determined 
according to 



K(v) = $>L(r)^W^(r) (16) 

a 



The vector d % a is a polarization vector of cx-th mode. In a crystal with cubic 
symmetry all modes can be classified as longitudinal or transverse ones, mak- 
ing the direction of the polarization vector independent of the orientation of 
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the elementary cell. Therefore, we can assume that local inhomogeneities, en- 
tering Eq.fll6|) through the r- dependence of the diagonalizing matrix c l aa (r), 
result in fluctuations of the effective charge Z a (r), leaving the directions 
of the polarization vectors without change. According to our basic assump- 
tion all triply degenerate modes contributing to the scattering have even 
parity, so the mean value of the effective charge is equal to zero. Correla- 
tion properties of the inhomogeneous effective charge are specified by means 
of the autocorrelation function K{r) = (Z cr (r)Z cr (0)) or its spectral density 
S(k) = J d 3 rK(r) exp(ikr). In what follows we assume that the correlation 
function, K(r), and spectral density S(k) are given as 



where (Z%) is an rms fluctuation of the effective charge, and k c is a correlation 
wave number, which is the inverse of a correlation radius and set the spatial 
scale of the inhomogeneities. We will also assume that the diagonalized 
dynamic matrix D a (r, r x ) is nonrandom, so that the mechanical properties 
of the system are not affected by the disorder. 

3.1 Vibrational response function 

Linear response functions T]^(Jl, q) and f]E{p<, q)can be found from the so- 
lutions of Eq. (|T5|) and Eq. (|T3|) with appropriate source terms on the right 
hand side of the equations. We will start by determining the mechanical 
response function rj^(Q, q). For this purpose we add an external force F to 
the right hand side of the equation Eq.(plJ), and our aim is to find an equa- 
tion for the averaged over the disorder amplitude, (W a (fl, K) . Carrying out 
the calculations with the use of a perturbation approach with respect to the 
rms fluctuation of the effective charge, cP, discussed in detail, for instance, 
in Ref. 1§, we can find with the accuracy to the first nontrivial correction 
(Bourret approximation) 



K(r) 
S(k) 




2 ■ 



(17) 



(fi 2 -4) (w a (n,K)- 

^-n 2 £<W a (n, K) / d 3 qS(K - q)4(K - qR(K - q)G'l(q) =F K . 

-oo » J 



(18) 
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We introduce here an initial dispersion law of the vibrations, e 2 K , describing 
triply degenerate (in the absence of interaction with electromagnetic field) 
mode. In the long wavelength limit the dispersion law assumes the form e 2 K = 
u!q + v 2 K 2 . Tensor G z J m (q) in Eq.([TJ) is the Green's function of Maxwell's 
equation fll3D in the absence of the coupling with the vibrations 

M 1 (x _ m \ i 1 VOi 

em[q) Q 2 -cV l j q 2 & q 2 ' 

High-frequency dielectric constant takes into account the electronic con- 
tribution into the polarization, and speed of light c in the previous equation 
is corrected for e^. Nondiagonal terms (with respect to the polarization in- 
dices a, a) in Eq.(pl|) can be shown to vanish after the integration over q, 
so we only have to evaluate r)Y for the longitudinal mode and rff for trans- 
verse modes. Evaluating the integral in Eq . (|T8|) one has to take into account 
that typical values of the correlation wave number k c determining the scale 
of the inhomogeneities can be of the order of magntidue betweenl0 6 cm -1 
and 10 8 cm _1 while the value of polariton wave numbers is of the order of 
10 3 cm _1 or 10 4 cm _1 . Another important observation is that the frequency 
u c = ck c corresponding to the correlation wave number is much greater than 
all other characteristic frequencies of the system. Therefore, one can neglect 
any effects of retardation, which are of the order of (0/ uj c ) 2 , and evaluating 
all the remaining integrals with the accuracy to (K/k c ) 2 , obtain longitudinal 
and transverse vibrational response functions in the following form 

*^ = iTH^Hai' (19) 

where 

4 = -o 2 + A 2 + (v 2 - l^)K 2 , 4 = .1 + (v 2 + i|)K 2 , (20) 

The parameter Tj„ takes into account initial damping of the vibrations, which 
is not related to the interaction with the field and is determined by intrinsic 
relaxation processes such as, for instance, anharmonism. This parameter can 
also include some effects of an additional scattering due to inhomogeneities 
that are intrinsic to the ion system, for example, spatial fluctuations of the 
dynamical matrix. Parameter A introduced above is an effective coupling 
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parameter defined as 

It can be seen from these expressions that the vibrational response func- 
tion does not contain effects of the disorder-induced interaction between vi- 
brations and transverse electromagnetic field. This happens because the 
effective relaxation parameter ck c , which is responsible for the scattering of 
the phonons into electromagnetic waves, is much greater than the effective 
coupling parameter A. Therefore, according to the general results obtained 
in Ref. 0, the coupling is ineffective and phonon properties are practically 
unaffected by the interaction. 

However, the interaction between phonons and a longitudinal macroscopic 
field modifies their properties. In particular, even zero-mean random effec- 
tive charge lifts the degeneracy between longitudinal and transverse phonon 
modes giving rise to the same LT-splitting that occurs in regular ideal crys- 
tals. It is more interesting that disorder also modifies the parameter v, which 
is responsible for the sign and strength of the spatial dispersion of the modes. 
The correction |p- has opposite signs for longitudinal and transverse modes, 
and its order of magnitude can be comparable to the initial parameter v. 
In ideal crystals the longitudinal electrostatic field is uniform and exerts an 
equal force upon ions at different sites, resulting in the uniform lift of the lon- 
gitudinal mode. In the presence of disorder, however, the field is not uniform 
anymore giving rise to certain corrections to the inter-ion potential. These 
corrections finally manifest themselves as modification of the parameter of 
the spatial dispersion v, and can be of the same order of magnitude as the 
initial potential causing substantial modification of the initial spectrum. 



3.2 Electromagnetic response function. 



In order to evaluate the electromagnetic response function r)£j(Q, K) one has 
to insert a source term into the Maxwell equation (|13|). According to the 
fundamentals of the linear response theory this source term should have the 
form of an external polarization ]nj]. Taking the source term into account 
and averaging the equations over the disorder in the same manner as above 
we finally arrive at the following equation for the average amplitude of the 
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electromagnetic field 
Air 



n 2 - c 2 k 2 



n 2 / d 3 q 



S(K - q) 



8 i:j + c 2 KiK 



1 3 



4ttQ 2 



(21) 



The electromagnetic response function rj?(Q, K) is usually determined as [JTT| 



(22) 



where f3ij(fl, K) is the inverse matrix for the left-hand side matrix of Eq.|2T|and 
is given as usual by the sum of the transverse and longitudinal parts 



PyfaK) = f3 T {5ij - KiKj) + (3 L KiKj, 



(23) 



where K is a unit vector in the direction of K.The integral with the spec- 
tral density given by Eq.([T7|) can be easily calculated to yield the following 
expression for (3 T 
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Q 2 - e\ - x 2 (VL) - 2itt(T in + T s 



{Q 2 - c 2 K 2 ) {Q 2 -e\- x 2 {Q) - 2iQ{T in + T s )) - A 2 Q 2 



(24) 



where 

x 2 (ft) = v 2 k 2 c -2vk c {l 



UJ, 



2\ 2 



1-7^ + 



r,(fi) 



OJ, 



2\ 2 



1 + 

Q 2 1 



Q 2 



AT 2 

in 

n 2 



4T 2 



-I 1/4 



sm 



- arctan 
2 



n 2 -u 2 



-I 1/4 



COS 



1 

- arctan 



2QT i 



n 2 



UOr 



(25) 



and we again introduced an initial phonon damping Ti n . The poles of the 
function (3t corresponds to the dispersion law of the transverse electromag- 
netic waves modified due to the random interaction with the phonons. The 
equation for the poles 

(SI 2 - c 2 K 2 ) (n 2 -e\- x 2 (fi) - 2ifi(T in + T s )) - A 2 ft 2 = (26) 
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resembles, as it was pointed out in Ref. |J, an equation for eigen frequencies 
of two waves propagating in a deterministic "effective" medium. The func- 
tions x 2 (f2) and r s (0) describe a renormalized dispersion law of phonons 
in this medium and their damping respectively. Actually, this equation de- 
scribes the interaction of the coherent (average) electromagnetic wave with 
non coherent (scattered ) ion vibrations, and its solutions give a dispersion 
law for the coherent component of the electromagnetic waves modified by the 
disordered coupling with the phonons. An analysis of the simplified version 
of the Eq.( p6|) carried out in Ref.|| [7|]) showed that this equation can yield 
a solution with two well-defined branches of " disorder- induced" polaritons if 
A > A cr , where A cr was estimated in Ref. JFj as A cr ~ v 2 k 2 /uo in the model 
without the initial relaxation T in . However, because of the small value of 
the parameter v of the phonon spatial dispersion, the initial phonon damp- 
ing cannot be actually neglected, and a more realistic expression for A cr can 
be shown to be A cr ~ a/ (T in /Q) (yk c + \fVtY in ) . Numerically the order of 
magnitude of A cr does not differ from the estimate presented in Ref.0, so 
the optimistic conclusion of that paper regarding the possibility of observing 
disorder-induced polaritons in real systems still holds. 

Two branches of the disorder-induced polaritons can manifest themselves 
as two maxima on the frequency dependence of the transverse component of 
electromagnetic response function r]fj(fl,K). Detailed study of the relation- 
ship between maxima of response functions and solutions of the dispersion 
equations was carried out in Ref. || in the case of the disorder-induced mag- 
netoelastic resonance. Applying the general results obtained in Ref.JJ to 
our situation we can conclude that two maxima appear on rjfj(Q, K) when 
a gap between two solutions of the dispersion equation (|26| ) at the crossing 
point becomes greater then A cr . Since A cr is the inversed decay time for the 
solutions at the resonance, the stated condition for two maxima to appear 
has a clear physical meaning - the frequency interval between two peaks must 
be greater than their width. 

Though qualitatively the conclusions drawn from the simplified model of 
Ref. |J remain valid for our more complicated situation, there is an impor- 
tant difference. The effective damping parameter in our case is frequency 
dependent. This fact leads to an asymmetrical form of the frequency de- 
pendence of the respective response function. From Eq.(|25|) one can see 
that parameter T s decreases with the increase of frequency. One can expect, 
therefore, that the low frequency maximum of the responce function will be 
lower and wider than the maximum corresponding to the upper polariton 
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branch. The plots presented in the subsequent section of the paper support 
this conclusion. 

The fact that vibrational and electromagnetic response functions demon- 
strate such a different behavior is typical for systems with disorder induced 
crossing resonances. The simple qualitative interpretation of this fact was 
given in Ref. ||, where the concept of two effective media was introduced. 
According to this concept, averaging over the disorder results in two differ- 
ent effective media for each of the interacting waves. The most important 
difference between the media is in their relaxation properties. The order of 
magnitude of the effective relaxation parameters is determined, among other 
things, by the slope of the initial dispersion curves and the correlation ra- 
dius of the disorder. The former in the case of electromagnetic waves is just 
speed of light and is much greater than the corresponding parameter v for 
phonons. As a result, a dispersion law for coherent phonons is determined by 
the interaction with electromagnetic waves propagating in strongly absorb- 
ing effective media; the modification of the initial dispersion law is negligibly 
small in this situation. As for coherent electromagnetic waves, their dis- 
persion law is modified due to the coupling with phonons in the effective 
media with rather weak absorption. The coupling is much more effective in 
this situation. Therefore, the initial dispersion curve splits at the resonance 
resulting in two maxima on the electromagnetic response function. 



4 Frequency dependence of the scattering cross 
section 

The response functions presented in the previous section can be used to 
study Raman scattering cross- section in our system. The results for the 
coherent part of the cross-section are obtained by a substitution of imaginary 
parts of Eq. (|T9| ) into Eq. (||D for the ionic contribution into the coherent cross- 
section, and the imaginary part of Eq. fl22|) into Eq. (^) for the electromagnetic 
contribution. As a result, the ionic contribution is 

jcoh _ 91,21 11 \ s s in iriQp 

j k m k p 2 (<W _ K m k p ^ 

I [W - el(K)} 2 + + [W _ 4 W ] 2 + AWT* 
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and electromagnetic part of the cross-section has the following form 



jcoh _ 2i 1 1 I s p s in i n 87rfiA 2 (r in + r 



2 («J mp - K m K^ Im/^(0, K) + K m i^Im/^(ft, #)] . (28) 



The Im/?^(f2, K) term in Eq.(p8|) represents the contribution due to scattering 
from transverse modes 

2fi 4 

lm(3r ~ [(n 2 - c 2 k 2 ) (n 2 - e 2 K ) - A 2 n 2 } 2 + 4n 2 (n 2 - c 2 K 2 ) 2 {v in + v s ) 2 ' 

(29) 

whereas the second term corresponds to the scattering due to the longitudinal 
mode 

lmp ' L = (fi 2 -4-A 2 ) 2 + 4fi 2 (r m + r s )^ ( 3 °) 

where we neglect the term x 2 (f2) defined in Eq.(^), which can be shown to 
be small. 

The polarization dependence of the coherent scattered intensity is mostly 
determined by the symmetry of the underlying structure, therefore we shall 
be focused upon a frequency profile of the intensity, which has more univer- 
sal significance. We would like to note, however, that the suggested average 
symmetry of the system results in vanishing of both coherent contributions 
when the polarization of the scattered wave coincides with the polarization 
of the incident wave. This feature of the considered structure can be used 
to separate the coherent contribution from the non coherent one. In what 
follows we restrict our consideration to the perpendicular polarization ge- 
ometry. In this case the angular dependent factors in Eqs. (|28| , p7|) take the 
form 1 — K 2 n for the transverse component , and K 2 n for the longitudinal 
one, where Ki n is a component of K in the direction of propagation of the 
incident wave, Ki n w sin(#/2), where 9 is an angle between wave vectors of 
the incident and scattered waves. 

Polariton effects manifest themselves in Raman scattering only at small 
scattering angles, corresponding to small wave numbers K. The main contri- 
bution into the coherent scattering in this situation comes from the transverse 
components, which depends upon the scattering angle as cos 2 (#/2), while the 
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contributions from the longitudinal components is small as sin 2 (#/2). The 
frequency profile of the spectrum is determined by the sum of two parts of the 
scattering cross-section: ionic contribution Eq. (|27| ) and electromagnetic con- 
tribution Eq. (|28|) . The former component of the cross-section does not con- 
tain any traces of polariton effects as was explained in the previous sections of 
the paper. Its frequency dependence has a simple form with one maximum, 
which is situated at the frequency of TO phonons. Due to weak phonon 
dispersion, position of this maximum almost independent of scattering angle 
6. Effects of disorder-induced polaritons can appear in the scattering spec- 
trum only due to the electromagnetic contribution to the cross-section. If the 
system under study allows for disordered polaritons to arise, the transverse 
response function Im/3^ has two peaks at the polariton frequencies, which are 
shifted with respect to the TO frequency. The resulting coherent intensity 
has, in this situation, three peaks: one in the center due to the ionic con- 
tribution, and two others due to the disorder- induced polaritons (see Fig.l, 
where electro-optical and mechanical coherent intensity are plotted together. 
A small bump on the electro-optical spectrum corresponds to a small lon- 
gitudinal contribution at LO frequency.) Positions and hight of polariton 
peaks significantly depend upon polariton wave number K, which in turn is 
determined by the scattering angle. These dependence can be obtained on 
the bases of the general analysis of susceptibilities presented in Ref. 0. At K 
much smaller than the resonance value, one can observe only one wide peak 
of Im/3r at LO frequency. The total coherent spectrum in this case will con- 
sist of two maxima corresponding to LO and TO phonons. Upon increasing 
K the original LO maximum will start shifting toward high frequencies while 
becoming narrower and taller. Simultaneously an additional maximum will 
emerge at frequencies below the frequency of TO phonons. The overall pic- 
ture in resonance region, hence, will contain three maxima with the central 
one corresponding to TO phonons, and two others resulting from disorder- 
induced polariton excitations. According to the general results of the theory 
of disorder-induced crossing resonance || the width of the polariton peaks al- 
lows for tracking the dependence between location of the peak and polariton 
wave number K, i.e. the dispersion law of the disorder-induced polaritons. 
If the crossing resonance does take place one would observe in this case the 
parts of two polariton branches described in Ref. 0, [jj . Increasing K beyond 
the resonance value will move the higher frequency maximum toward even 
higher frequency while making it sharper. The low frequency maximum at 
the same time will be approaching TO frequency and broadening. Finally 
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the upper polariton branch will become unvisible, and the lower one will 
be indistinguishable from TO phonons. At the same time the contribution 
from longitudinal cross-section will rise, therefore one will again observe the 
spectrum with two maxima at LO and TO frequencies. 

If the effective coupling parameter A is too small for the disorder-induced 
crossing resonance to occur, then all one could observe is a two maximum 
spectrum whose shape only weakly depends upon the scattering angle. 

Now let us consider non coherent scattering. Substituting Eq.(|19D for the 
vibrational response function in Eq . (|10|) for the non coherent cross-section, 
and evaluating the integral with the spectral density given by Eq.([H]), one 
arrives at the following expression for I^ oh 
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and v\ L = (v 2 ± |A 2 //c 2 ) are the phonon dispersion parameters modified 
due to the interaction with the electromagnetic field (see Eq. (POf) ) . we have 
neglected here if- dependence of the phonon dispersion laws. Eq.([H]) has a 
structure similar to the expression for the coherent vibrational intensity ex- 
cept for its independence from the scattering angle. Therefore, transverse and 
longitudinal components have equal weights. The positions and the widths of 
the both peaks, however, differ from the coherent scattering. Fluctuations of 
the mechano-optical coefficients result in a shift of the transverse and longi- 
tudinal frequencies by v\ L k 2 c respectively, and cause an increase in the width 
given by the second term in Eq . (p2|) . This expression is similar to Eq . (p5|) , 
which determines the widths of the maxima of the coherent electromagnetic 
contribution, with the obvious replacement of the parameters. 
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The analysis of the non-coherent electromagnetic scattering is somewhat 
more cumbersome because of the more complicated polarization dependence. 
We would like to point out two simplifying details. First, a contribution from 
the transverse part of the electromagnetic response function has the order 
of (il/ck c ) 3 and can be neglected. Second, we can neglect independence 
in the corresponding integrals. As a result, the expression for non coherent 
electromagnetic part of the scattering intensity takes the following form 
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where G\ and G2 are certain combinations of the original coefficients G. One 
can see that the non coherent electromagnetic contribution is identical to the 
coherent longitudinal contribution. 

A sample plot for the total scattering intensity is shown in Fig.2. In gen- 
erating plots in Fig.l and Fig.2 we choose parameters typical for the systems 
under consideration, the most important of which are the coupling parame- 
ter, internal phonon damping rate, electro-optical and mechano-optical coef- 
ficients. The coupling parameter A was chosen to be an order of magnitude 
smaller than the resonance frequency, and phonon damping was assumed to 
be two orders of magnitudes smaller than the frequency. The relative value 
of the electro-optical and mechano-optical parameters was chosen to make 
respective contributions into the scattering intensity be of the same order 
of magnitude. One can see that the total scattering retain the three-peak 
feature of the coherent contribution. The peak at the center is formed as a 
result of superposition of the mechanical coherent cross-section and all the 
non coherent contributions. Obviously the width of the non coherent max- 
ima for values of parameters chosen for the plot appears to be to large to 
reveal different types of non coherent contributions, which formally should 
have appeared at frequencies shifted with respect to coherent TO and LO 
frequencies. Since the values of the most parameters required to make more 
explicit predictions of the shape of the spectra are not known we cannot 
discuss details of the spectra. The figures presented demonstrate neverthe- 
less that under favorable circumstances non coherent contribution would not 
mask effects due to disordered polaritons in Raman scattering intensity. We 
can conclude, therefore, that one could hope to observe the polariton- related 
effects in materials whose average structure forbid polariton formation by 
studing the spectra for different scattering angles as explained above. 
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5 Conclusion 



In this paper we develop a theory of Raman scattering of light due to 
disorder-induced polaritons. We consider a material that on average retains 
a diamond-like cubic structure. This structure has only one triply degener- 
ate optic phonon mode of even parity, which is dipole-inactive in an ideal 
structure. Local random distortion of the structure, however, gives rise to 
a random effective charge of the mode and respectively to an effective ran- 
dom coupling with electromagnetic field. The distortion can be induced by 



different means such as application of pressure |12[], using thin films fT3|| . 
or amorphization. Among materials with the structure of diamond, for ex- 
ample, Si and Ge exist in amorphous versions. Raman scattering in the 
systems considered in the paper has an interesting feature, which does not 
show up in the case of ideal structures. In crystalline materials with center 
of symmetry there exists an "exclusion rule" for dipole infra-red activity and 
Raman activity. Dipole-active modes, which form polariton states, do not 
participate in the first-order Raman scattering. Therefore, in ideal crystals 
the problem of polariton Raman scattering arises only in structures without 
a center of symmetry. In our situation we deal with a mode which is, on av- 
erage, dipole inactive, and is, therefore, Raman active. At the same time this 
mode is coupled to electromagnetic field due to fluctuating effective charge, 
and induces random electromagnetic field. This field has zero average am- 
plitude, but its intensity is not zero, and it makes sense to take into account 
the modulation of the electron polarizability due to this field. Components 
of the respective tensor of the susceptibility derivatives are, of course, ran- 
dom functions and contribute to the non coherent scattering. The average 
symmetry of the structure does not forbid, however, this tensor from having 
non- zero average value. If this is the case, then one can observe coherent elec- 
tromagnetic component of the scattering cross-section from disorder-induced 
polaritons. Assuming that the average electro-optical coefficients in our situ- 
ation do not differ in magnitude from more standard situations we show that 
polariton contribution into Raman scattering can indeed be observed on the 
non-coherent background. 

Raman spectra of amorphous materials have been studied for many years 
(see, for example, the latest paper in Ref|l4| and references therein). Their 
interpretation, however, is rather complicated, therefore, it is difficult to de- 
termine if the disorder-induced polaritons have been already observed. An 
accurate identification of disorder-induced polaritons in the spectra requires 
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special experiments. We considered the behavior of the Raman cross-section 
under the change of the scattering angle, and determined which features of 
the spectra are specific for disorder-induced polaritons. On the basis of our 
analysis we think that comparison of spectra measured at different scatter- 
ing angles could reveal the effects caused by disorder-induced polaritons. We 
suggest that crystals with locally distorted structures rather than amorphous 
materials are the best candidates for disorder-induced polariton experiments. 
Additional means to interpret experimental Raman spectra are provided by 
the possibility to compare spectra of disordered materials with their crys- 
talline counterparts. In this case one has to take into account that our model 
does not explicitly describe the shifting and broadening of the spectra due 
to inhomogeneities intrinsic to the phonon subsystem. We incorporated the 
broading phenomenologically by means of the initial damping parameter Ti n , 
and can also assume that initial phonon TO frequency is already corrected 
for the shift due to the intrinsic phonon scattering. Such consideration is 
valid if one assumes that intrinsic phonon parameters of the system (density 
and force constants) fluctuate independently of the phonon-photon coupling 
parameter. This is a legitimate assumption since a correlator between, for in- 
stance, the force constants and the coupling parameters should be the tensor 
of rank three. Such a tensor must vanish in a system with the average struc- 
ture displaying a center of symmetry. Therefore, one can say that the central 
peak of the predicted three-peak spectrum corresponds to the broaden and 
shifted peak, which would correspond to TO phonon peak in the abscence of 
the polariton effects. 

Among other results we would like to mention the modification of phonon 
dispersion laws due to the fluctuating longitudinal electromagnetic field. Our 
calculations show that in random systems the longitudinal field modifies 
both the fundamental frequency of phonons and the slope of their dispersion 
curves. The last effect means that fluctuating component of longitudinal 
field contributes into interactions between atoms from different elementary 
cells. It is interesting, that this contribution is found in transverse as well as 
in longitudinal modes of the vibrations, but with opposite signs. The fluc- 
tuating field weakens longitudinal components of the interatomic forces, and 
enhances the transverse ones. 
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Figure Caption 



Fig.l Double peacked line shows electro-optical contribution to coherent scat- 
tering intensity. The single peak line is the coherent mechanical contribu- 
tion. The frequency on the plot is normalized by the resonance frequency. 
Fig. 2 Total scattering intensity with three peaks. The left one presents the 
contribution from the lower polariton branch, the right one from the up- 
per polariton branch. The peak in the middle is the resulting contribution 
from TO phonons and non-coherent scattering. Upon increase of the non- 
coherent scattering the left peak will dissapear first, and yet the remaininmg 
double-peaked structure will present a clear indication on the existence of 
the disordered polartions. 
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